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ABSTRACT 

We present detailed calculations of the magnification distribution, including both weak 
and strong lensing, using very recent solutions of the Dyer-Roeder (1973) equation for light 
propagation in a inhomogeneous universe with a cosmological constant and up-to-date mod- 
els for the evolving cosmological distribution of dark matter halos. We address in particular 
the effect on the magnification distribution of a non-zero cosmological constant, of different 
density profiles of lenses, and of finite sizes of lensed sources, three important issues not yet 
fully settled. 

We show that, if dark matter fluctuations are normalized to the local cluster abundance, 
in the presence of a cosmological constant the optical depth for lensing decreases compared 
to the case of an Einstein-de Sitter universe, because halos in the relevant mass range are 
less abundant over a large redshift interval. We also discuss the differences in the magnifica- 
tion probability distributions produced by Navarro, Frenk & White (NFW) and by Singular 
Isothermal Sphere (SIS) density profiles of lenses. We find that the NFW lens is more efficient 
for moderate magnifications (2 ,$ A & 4), and less efficient for larger magnifications. More- 
over, we discuss quantitatively the maximum magnification, A max , that can be achieved in the 
case of extended sources (galaxies) with realistic luminosity profiles, taking into account the 
possible ellipticity of the lens potential. We find that plausible values of A max are in the range 
10-30. 

Finally, we apply our results to a class of sources following the luminosity evolution typical 
for a unified scheme of QSO formation. We find that the lensed source counts at 850/jm can 
be larger than the unlensed ones by several orders of magnitude at flux densities lOOmJy. 



1 INTRODUCTION 

Gravitational lensing is well known to be a powerful tool to probe 
the overall geometry of the universe at z ~ 6 and the values of cos- 
mological parameters (curvature, vacuum energy density described 
either by a cosmological constant or by a dynamical quantity like 
quintessence, Hubble constant: Bartelmann et al. 1997; Falco et 
al. 1998; Cooray 1999; Huterer & Cooray 1999; Macias-Perez et 
al. 2000; Bhatia 2000; Helbig 2000), the evolution of large scale 
structure (Rix et al. 1994; Mao & Kochanek 1994; Bacon et al. 
2000), the masses and the density profiles of dark halos of galaxies 
and galaxy clusters (Narayan 1998 and references therein; Mellier 
1999; Clowe et al. 2000). 

Another important effect of gravitational lensing is the modi- 
fication of the observed luminosity functions of distant sources and 
of their counts, due to the redshift dependent magnification bias 
(Peacock 1982; Narayan 1989; Schneider 1992). 

In view of all the above applications, it is essential to perform 
the most accurate modeling of all the ingredients entering into the 
gravitational lensing process. Quantitative predictions depend on 
the lens population as a function of mass and redshift, on the lens 



profiles, on the the luminosity function and size of the sources, and 
on a correct treatment of the distance-redshift relations, taking into 
account the effect of density inhomogeneities on the propagation of 
light rays. 

Gravitational lensing by intervening mass concentrations is a 
statistical process. It can be described by a probability distribution 
F(A,z s ) of the magnification factor A, which obviously depends on 
the source redshift z s . It is generally difficult to model F(A,z s ) in 
detail analytically. Some authors estimated it with A'-body simula- 
tions (see e.g. Rauch 1991; Pei 1993a,b). Nonetheless, the main 
features of F(A,z s ) can be captured using generic properties of 
gravitational lensing. Strong lensing occurs near caustic curves, 
which cause F(A,z s ) to decrease oc A~ 3 for A> 1, irrespective 
of the lens model (Peacock 1982; Turner et al. 1984). 

In this paper we present detailed calculations of the magnifi- 
cation distribution, including both weak and strong lensing, using 
very recent solutions of the Dyer-Roeder (1973) equation for light 
propagation in a inhomogeneous universe with a cosmological con- 
stant (Demianski et al. 2000) to compute the relationships between 
distances and redshift as well as up-to-date models for the evolving 
cosmological distribution of dark matter halos (Sheth & Tormen 
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1999; Sheth et al. 1999) that can act as lenses. We address in par- 
ticular the effect on the magnification distribution of a non-zero 
cosmological constant, of different density profiles of lenses, of a 
finite size and of luminosity profiles of lensed sources, three im- 
portant issues not yet fully settled. The magnification probability 
distribution has been applied to the counts at 850/jm: the luminos- 
ity function of the source population used here has been inferred 
by the evolution of the QSOs luminosity function and from the link 
between the QSOs and the hosting spheroids (Granato et al. 2000). 
Such luminosity function turns out to fit the SCUBA data. 

The plan of this paper is as follows. In Sect. | and | we de- 
scribe our assumptions on the distribution of lenses in mass and 
redshift, and on their density profiles, respectively. In Sect. ^[ we 
review the basic definitions of distance and flux magnification in 
inhomogeneous universes. In Sect. ^, we introduce the probability 
for a source at given redshift to be magnified by a certain amount, 
and quantify the magnification bias on source counts, for which 
a model is described in Sect. ^ In Sect. ^ we face the problem 
of determining the maximum magnification allowed for extended 
sources. We summarize and discuss our results in Sect. || 

Throughout this paper, Q.Q m and £2oa denote the present-day 
density parameters for the non-relativistic matter and for the cos- 
mological components, respectively, neglecting the radiation en- 
ergy density. The Hubble constant is Hq = lOO/zkms -1 Mpc~', 

= 1 — flom — &oa ls the curvature term, and fl denotes the to- 
tal density parameter. Our fiducial model has flom = 0.3, £2qa = 0.7 
and h = 0.65. Furthermore, qg = —da/d 2 is the deceleration param- 
eter, and we concentrate on spatially-flat models (O.o^ = 0). 



2 THE EVOLVING HALO MASS FUNCTION 

We assume that the lens population consists of collapsed dark mat- 
ter halos with an epoch-dependent mass function described by the 
Sheth & Tormen (1999) function, which reproduces fairly accu- 
rately the results of extensive numerical simulations over more than 
four orders of magnitude in mass, for a wide range of Cold Dark 
Matter cosmologies (Jenkins et al. 2000). This function improves 
considerably on the familiar Press-Schechter (1974) model which 
overestimates the abundance of "typical" (M*) halos and underes- 
timates that of massive systems. The comoving number density of 
halos with mass M at redshift z is then 
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The best-fit values of the parameters are a = 0.707, p = 0.3, and 
A ~ 0.3222 (Sheth & Tormen 1999; Sheth, Mo & Tormen 1999). 
The Press-Schechter mass function is recovered for a = 1, p = 
and A = 0.5. 

In Eq. ([!]), po is the mean mass density at a reference epoch 
fo, which we assume to be the present time, and a 2 is the variance 
of linear density fluctuations at the present epoch smoothed with a 
spherical top-hat filter Wr (k) enclosing a mass M. The variance a 2 
is related to the power spectrum P(k) by: 

471 r &kk 2 wi(k)P(k) , (2) 
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A useful expression for a(M) and its derivative dc/dM was ob- 
tained by Kitayama & Suto (1996), 

o°c (l+2.208m rf -0.7668ra M +0.7949m 3£/ r 2/(,;W) , (4) 



r = Q.o m hexp[—Q.Qi > ({ + v2h/Q.o m )] being the shape parameter 
(Bardeen et al. 1986) and we have adopted a present baryon density 
Sloth 2 = 0.03. For a ACDM model with Cl 0A = 0.7 and il m = 0.3, 
the COBE normalization is Og = c(R = 8/i~' Mpc) = 0.925 (Bunn 
& White 1997). 

Also in Eq. ([!]), S 2 (z) is the linear density contrast of an object 
virializing at z, evaluated at the present epoch. It can be estimated 
using the spherical collapse model (e.g. Peebles 1980; Lahav et 
al. 1991; Lacey & Cole 1993; Nakamura 1996; Eke, Cole & Frenk 
1996; Lokas & Hoffman 2000) yielding: 



5c (z) 



8 c D + (z = Q) 
D+(z) 



(5) 



where D + (z) is the linear growth factor of density fluctuations: 



n+zr'^n 



x [^ /7 -fl A + (l+n m /2)(l+^ A /70)]- 1 , (6) 

in the approximation by Carroll et al. (1992), where f2 m = ilom(l + 

z) 3 /E 2 (z),Cl A = noA/E 2 (z),tmd 

E(z) = H(z)/H = [fl 0m ( 1 + z) 3 +n 05l (1 + z) 2 + £2oa] 1/2 . (7) 

8 C [Eq. (|J)] is the value of the critical overdensity evaluated at viri- 
alization. An approximate expression valid for any cosmological 
model was reported by Nakamura (1996) and Kitayama & Suto 
(1996): 

3(127i) 2 / 3 

8c « K 2Q > (1 +0.0123 log 10 n m ). (8) 



3 GRAVITATIONAL LENSING AND HALO MODELS 

The ray-tracing equation relates the position of a source to the im- 
pact parameter in the lens plane of a light ray connecting source 
and observer. The light ray passing the lens at an impact parameter 
t, is bent by an angle 6l(Q. The source position f) and the impact 
parameter £ in the lens plane are related through 



*i = 7rC-Ai S &(C) , 

At 



(9) 



where Dd.s.ds are the angular-diameter distances between observer 
and lens, observer and source, and lens and source, respectively. 
We will specify these distances later in the context of a "clumpy 
universe". The deflection angle a is the sum of the deflections due 
to all mass elements of the lens projected on the lens plane, 



4G 



Z(C')d 2 C, 



(10) 



where Z is the surface mass density of the lens. 

It is convenient to introduce the dimensionless impact parame- 
ter x= C/Co> where i^o is an arbitrary length scale. The correspond- 
ing value projected in the source plane is T|o = i^o^s/Dd an d the 
source position is y = fj/r|o. 

The dimensionless surface mass density k(x), also called con- 
vergence, is 
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k(x) = — — with Z cr = 

V ' £cr 4jiG D d D ds 

The dimensionless ray-tracing equation is then, 

y = x — 6c(x) , 

with the scaled deflection angle 

a(x) = ^£^a(c, x) = - f dx'x'K(x') . 

The axial symmetry of the lenses considered here allows us to 
rewrite the lens equation in terms of the dimensionless "mass" m(x) 
enclosed within the radius x, 



(12) 



(13) 



m(x) 



(14) 



with 



mix) =2 \ dx'x'K(x') . (15) 
Jo 

The Jacobian matrix H (x) = dy/dx for the lens mapping [Eq. (|l4|)] 
determines the magnification factor for each image i, 

« = . • (16) 

det Jl{xi) 

The sign of /j; reflects the parity of the image with respect to the 
source. The total source magnification is the sum of the absolute 
values of the magnification factors of all images, 



(17) 



For axially symmetric lenses, the Jacobian determinant of the lens 
mapping can also be written 

a(x) s 

X 

d / .'.'/ n 

dx v x , 



det J? 
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where yis the shear (cf. Schneider, Ehlers & Falco 1992). Critical 
curves are located at the zeros of det;3 ; their images in the source 
plane under the mapping described by Eq. |T^ ) are the caustics. 

A simple model for the mass profile of a lens (cluster or 
galaxy) is the singular isothermal sphere (SIS; e.g. Binney & 
Tremaine 1987): 

P «=2^' (19) 
where o v is the line-of-sight velocity dispersion. In this case, the 
deflection angle is independent of the impact parameter (cf. Schnei- 
der et al. 1992; Narayan & Bartelmann 1997). One of the two crit- 
ical curves degenerates to a point. Any given source has either one 
or two images. Two images appear only if the source lies inside the 
Einstein ring. The Einstein angle corresponding to a SIS is 

,»ds 



8f. = a - 



with 



:4ji^t = 1.4" 
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220kms~ 



(20) 
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In order to evaluate the mass of a halo with velocity dispersion a v 
from Eq. (|l9[), we truncate the sphere at an effective radius com- 
puted following Lahav et al. (1991). 

If every halo virializes to form a singular isothermal sphere, 
mass conservation implies that the velocity dispersion is related to 



the mass and to the virialization redshift by (see, e.g. Porciani & 
Madau 2000; Kaiser 1986): 



T 1/6 



,1/2 



(22) 



where Q. m is defined in the previous section, ro = (3M/4jipo)'' 3 , 
and A(z) is the mean density of the virialized halo in units of the 
critical density at that redshift. A useful expression for the depen- 
dence of o v on redshift is given in Bryan & Norman (1998), 



o,=Jf 1 /3[fl 2 (z)A(z)G 2 /l6] 1/6 



(23) 



Bryan & Norman (1998) also give a fitting formula for A(z) for a 
flat universe with a cosmological constant and for a universe with 
^Ost 7^ an d = 0. For a flat universe, (£2o^ = 0) their hydro- 
dynamical simulations yield: 

A(z) = 18tt 2 + 82x - 39x 2 , (24) 



where x = £l m — 1 . 

The SIS model is useful because it allows to work out an- 
alytically the basic lensing properties. On the other hand, high- 
resolution A/-body simulations (Navarro, Frenk & White 1997) 
showed that in hierarchically clustering universes, virialized dark 
matter halos have a universal density profile (referred to as NFW), 
which is shallower than isothermal near the center and steeper in 
the outer regions: 



p(x) 



PcritONFW 

x(l+x) 2 



(25) 



where x = r/r s , r s being a scale radius, 8nfw is the characteristic 
density contrast of the halo, and p cr ; t is the critical density at the 
epoch of the halo virialization. This formula was found to accu- 
rately describe the equilibrium density profiles of dark matter halos 
over a broad range of masses (3 x 10 11 < M m /M @ < 3 x 10 15 ), 
irrespective of the cosmological parameters and the initial density 
fluctuation spectrum. 

Still higher resolution A/-body simulations indicate a steeper 
central cusp than that of the NFW profile: p(x) <* [# (1+x) 1 
(Moore et al. 1999; Ghigna et al. 2000). Since the slope of the cen- 
tral density profile in this case [p(r) o= r~ 1,5 ] is intermediate be- 
tween those of the NFW [p(r) °c r~'] and of the SIS [p(r) °c r~ 2 ], 
the two cases considered here (NFW and SIS) will bracket it. 

We parameterize halos by their mass M200 (defined as the 
mass within the virial radius 7-200, the radius of a sphere of mean 
interior density 200p cr i t ; see Navarro et al. 1997), since this is 
the mass whose distribution is given by Eq. (jl|). Unless otherwise 
specified, we abbreviate M200 by M. The halo concentration is 
c = ^20o/ r s- It is related to the density parameter S^pw by 



200 



3 pn(l + c)-c/(l+c)] 



(26) 



The virial radius of a halo at redshift z depends on the halo mass as 



^200 : 



1.63 x 10" 
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/r'kpc. (27) 



Alternatively, the halo can be characterized by its circular velocity, 
V200 = f^)' /2 



^200 

h~ l kpc 
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(28) 
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The scale radius r s depends on the halo mass. The halo concen- 
tration increases with decreasing halo mass (e.g. Fig. 6 of Navarro 
et al. 1997). Less massive halos are therefore more concentrated. 
For a given halo mass, Eqs. ( p7| ) and ( po| ) completely specify the 
density profile [Eq. (|2^)]. 

The lens equations for the NFW profile are given by Bartel- 
mann (1996) and Maoz et al. (1997). The surface mass-density is 



£ M = x 2 _ ! /to- 



with 



/to = 



1- 

1- 




. arctan 



=i=arctanh 
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The dimensionless lens mass is 

, \ 4p cl -i t 8 NFW r s 
m{x) = g(x) , 



with 



g(x)=ln- + l-/(x) 



(29) 



(30) 



(3D 



(32) 



While the NFW profile [Eq. (J25[)] has two critical curves (Bartel- 
mann 1996), the SIS profile has only one. A SIS lens has either one 
or two images, an NFW lens either one or three. 



4 THE MAGNIFICATION BIAS 

In this section, we follow the approach by Schneider et al. (1992) 
and Schneider (1987a,b). We call /j the magnification of fluxes in 
an inhomogeneous universe, and relate it to the magnification A in 
a homogeneous universe. 

The effect of lensing on flux-limited source counts is quanti- 
fied by the magnification bias (Turner et al. 1984). If the surface 
density of galaxies with flux greater than S v is N(S V ), flux magni- 
fication alters the counts to N'(S V ) = j" -1 -/V(Sv/J~' ), and the mag- 
nification bias is given by (Narayan 1989): 



(33) 



The factor arises because solid angles are magnified, hence 
source counts are diluted. Eq. (B3J) implies that for power-law 
counts (N oc Sy a ), q(/j,S v ) = jU , so that the magnification bias 
increases as the counts steepen. 

Since the reference to "unlensed" counts has given rise to con- 
fusion in the literature, it is important to specify our model for 
light propagation in the universe when talking about the number 
of sources seen "in absence of lensing". 

In Eq. (j33|), N is defined as the number count of sources that do 
not appear behind lenses. Because of energy conservation, the flux 
of these objects must be reduced compared to a homogeneous uni- 
verse with the same average density. In other words, since lensing 
is caused by inhomogeneities, we have to deal with an inhomoge- 
neous universe. While sources observed through lenses are magni- 
fied by a factor /j + > 1 compared to a homogeneous universe, un- 
lensed sources must be demagnified by /j~ < 1 . When comparing 
lensed to unlensed sources, the effective magnification is 



(34) 



(e.g. Schneider 1984, 1987a,b). This is the "empty beam approach" 
to light propagation (Dyer & Roeder 1973; Ehlers & Schneider 
1986), in which light cones are devoid of clumped matter. The frac- 
tion of uniformly distributed matter is denoted by a s , the smooth- 
ness parameter. 

In this on-average homogeneous and isotropic universe, the 
average flux (5) from a source population at redshift z with lumi- 
nosity L, must be equal to the corresponding flux Sfl mat would be 
observed in a strictly uniform Friedman-Lemaitre universe, 



(S) = Sfl 



LK(L,z) 



(35) 



where Z)l (z) is the luminosity distance in the uniform universe, and 
K(L,z) is the K-correction. In general, it is not possible to compare 
mean fluxes of sources in these two different spacetimes. In par- 
ticular, the notion of distance has no unique meaning in a clumpy 
universe, as it depends on both redshift and direction. Due to the 
corrugated structure of the gravitational field, the propagation of 
light in the inhomogeneous universe is a statistical problem. We 
can compare fluxes as in Eq. ( p5[ ) only because we assume that the 
global geometry of the clumpy universe equals that of the homo- 
geneous universe. Equation (B5J) implies that the area of a surface 
of constant redshift on the future light cone of a source is equal 
to that of the corresponding wavefront in the Friedman-Lemaitre 
model. In order to use redshift-distance relations in an inhomoge- 
neous, partly clumpy universe, which is uniform on average, it is 
then necessary to introduce a meaningful operational definition of 
distance. Since light propagation in a clumpy universe depends on 
the clumps in and near the light beam, Dyer & Roeder (1973) con- 
sidered the limiting case of a light bundle that propagates far away 
from all clumps and is thus unaffected by gravitational shear. This 
is the "empty cone" limit. 

The angular-diameter distance is then computed by replacing 
the density parameter £l m with the "homogeneous" fraction a s fi m . 
The resulting "Dyer-Roeder" distance is larger than the angular- 
diameter distance in a homogeneous universe. For a spatially-flat 
universe with deceleration parameter qo = 0.5, the Dyer-Roeder 
distance between redshifts z\ and Z2 > Z\ is 



(l +Z2 )(P-5)/4 (l +zl )(0-5)/4 



(l+zi)(P+ 5 )/ 4 



(l+za 



)(P+5)/4 



(36) 



where p = (1 + 24£2 G ) 1//2 and Sl c = (1 - a s )fl m is the density 
parameter in compact objects (P = oc s = 1 for the homogeneous 
universe). Hereafter, we abbreviate r(0,z, P) = r(z, P) and r(z, 1) = 
r\ (z) ■ The dimensional angular diameter distance in a homoge- 
neous universe is 



Di(z) 



n (z 



(37) 



The problem is more complicated when the cosmological constant 
is positive. The Dyer-Roeder equation in a FRW universe with non- 
zero cosmological constant has been solved exactly by Kantowski 
(1998), Kantowski & Kao (2000) and, following a different ap- 
proach, by Demianski et al. (2000). The general solution can be 
expressed in terms of the hypergeometric functions f s ± introduced 
by Demianski et al. (2000): 



r(z,P)=Ai 



(l+z)- p/4 /* 



A 2 T< . WTa /, , _^ i (38) 



(1+Z)V4 (1 + Z )3 Z (l+ Z )5/4 (l+ z )3 



where the constants A 12 are determined by the boundary conditions 



© 0000 RAS, MNRAS 000, 000-000 



Gravitational lensing 5 



r( Z ,(3)| z=0 = 



and — 

dz 



1 



(39) 
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This can be used to generalize the two-point distance of Eq. 
which, for a non-zero cosmological constant, reads: 



Kzi>Z2)P) 



Kzi,P)(l- 



s: 



dz' 



fz') 3 ^0m4 



■OA 



1/2 



r2(z',p)(l+ Z ') 



A2 



(40) 



with r(z, (3) given by Eq. pq). 

For practical applications, these exact solutions are rather 
cumbersome to use. Demianski et al. (2000) give approximate so- 
lutions for r(z, P) and r(t\ ,Z2 , P), which we use in this paper in the 
form kindly provided by R. de Ritis (private communication) for 
Oqa = 0.7 and <x s = 0.9, implying P = 1.84, and Q G = 0.1. We ig- 
nore a possible time-dependence of the smoothness parameter. The 
effect of different choices of oc s on the Dyer-Roeder distance with 
positive cosmological constant is shown in Demianski et al. (2000). 

Having introduced the Dyer-Roeder distance, we can inter- 
pret the magnification in Eq. (B4J) as the ratio between the flux S 
actually received from a source at redshift z, and the flux Sempty 
that would be received if the same source was observed through an 
empty cone: 



= 5 



J empty 



4iiP£(z,P) 
LK(L,z) 



> 1 



where 



^teP) = flo(l- 



K«.P) 



(41) 



(42) 



is the luminosity distance in an empty cone. 

Since the clumpy universe is homogeneous on-average, 
Dl(z, P) = Dl(z, 1). Using Eq. @, Eq. @ reads 

2 

> 1 . (43) 



<A<> = ( 



rfeP) 



n (z) 



(fi) can be interpreted as the magnification of an average light beam 
in a smooth universe relative to that for an empty cone in a clumpy 
universe: (p.) = Sfl /Sempty 

For go = 1/2, (p) becomes 



(fi) = P~ 



sinh 2 [0.25pin(l+z)] 
sinh 2 [0.251n(l+z)] 



(44) 



(Pei 1993b), while for the more general case including the cosmo- 
logical constant, (//} has to be derived from Eq. 

Let p(p,z)dp denote the probability for a source at redshift z 
to be magnified by a factor of p within dp. The normalization and 
flux conservation conditions require: 

dpp{p,z) = 1 , / dpp{p,z)p = (p) . (45) 

When dealing with source counts on cosmological scales, however, 
it is often convenient to refer magnifications to the homogeneous 
universe rather than to the "demagnified" background. We call A 
the magnification relative to the homogeneous universe: 



S 
Sfl 



J empty 

Sfl 



(46) 



where p = S/S em pty is the magnification of a source seen through 
an empty beam, as in Eq. (plj). Eqs. (^5j) and (^) imply: 



A = 



(f) 



(47) 



This is the definition of magnification used, for example, by Blain 
(1996) and Peacock (1982). It is trivial to show that for a given 
source redshift z, p(p,z)dp = p(A,z)dA, and the normalization and 
flux conservation conditions become 



/: 



dA p(A,z) = 1 



dA p(A,z)A = 1 



(48) 



Here, A mm is the minimum magnification relative to the FRW uni- 
verse, that can be produced by compact lenses, i.e. A m ; n = {/./) ~~ 1 . 
Since (A) = 1, the sources which are demagnified with respect to 
the homogeneous universe have A < 1. 



5 THE MAGNIFICATION DISTRIBUTION 

The number of images produced by a lens, their angular separation 
and their magnification depend on their relative alignment as seen 
by the observer. For a fixed geometry of the lens system, one can 
ask where the source must lie for its images to have a total magni- 
fication larger than a given value p. The area of the resulting region 
in the source plane is the cross section a(p,Zd,Zs,%), which obvi- 
ously depends on lens and source redshifts, Zd and z s , on a set of 
parameters % describing the lens model, and on the magnification 
itself. As both SIS and NFW halos are completely characterized 
by their mass, % = M here. The cross section quantifies the effi- 
ciency of the individual lens on a source. It is generally found nu- 
merically, by solving the lens equation and finding the area in the 
source plane where a source must lie in order to produce images 
with a total magnification larger than /j to t- The cross section of a 
singular isothermal sphere for a total magnification /j t ot > M is 
4jj<3[2 q1 

o(p tot >p) = j-4 for p>2. (49) 

We need to compute the total magnification cross section of an 
ensemble of lenses distributed according to the mass function of 
Eq. (|l|). We denote by n c (z,M) the comoving lens number density, 
and n(z,M) = (1 +z) 3 w c (z,M). As long as the cross sections of 
individual lenses do not overlap, each light bundle from a source 
encounters only one lens, and the total cross section is the integral 
of the individual cross sections over the redshift and mass distribu- 
tions of lenses: 

3 



OtotC",Zs) 
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(50) 
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c(p,z,z s ,M)n c (z,M)(l 



^2 „2 



v/OOmO +z) 3 + fl A 
where the proper volume of a spherical shell of width dz at redshift 
z is 
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For a spatially-flat universe 
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dz. 



'prop 



(51) 



(52) 



d z tto(l+z)\AWl+z) 3 +CioA 
The probability for a source at redshift z s to be lensed with magni- 
fication > p is obtained by dividing o tot by the area of the source 
sphere 

r?(z)(l+z) 2 



i 



dz- 



Ho rf(z s ) JO v/^0m(l+z) 3 +fi0A 

J dMa(p,z,Zf,,M)"c{z.M) . 



(53) 
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The requirement of non-overlapping cross sections restricts the va- 
lidity of Eq. (J5|> to P <C 1, i.e. the total cross section must be much 
smaller than the area of the source sphere. The net effect of gravi- 
tational lensing on the distribution of flux densities expected from 
a population of distant sources can be described by the probability 
distribution of magnifications, p(p,z) introduced in Sect. Q While 
P< 1, the differential probability is p(p,z) = —dP(p,z)/dp; in par- 
ticular, p{A,z) = p(p,z)(p). 

The differential probability decreases as p~^ for p 2> 1, hence 
the high magnification tail in terms of A can be written as p(A,z) °^ 
a(z)A~ 3 . On the other hand, Eq. ( |53| ) breaks down for small 
magnifications, where multiple lensing events become important 
and cross sections overlap (in fact, the probability for many low- 
magnification lensing events along the line of sight to a source is 
rather large, while a single interaction producing high magnifica- 
tions is a relatively rare event). In particular, this implies that there 
is a critical magnification A cut below which multiple lensing be- 
comes important, resulting in low-magnification events (weak lens- 
ing regime). Vietri & Ostriker (1983) analytically described multi- 
ple lensing for the SIS model. 

Based on general considerations (see Bartelmann & Schnei- 
der 2001 for a review), one expects a skewed magnification dis- 
tribution with a weak lensing peak near A = 1, turning into the 
high-magnification tail °= A~ 3 at A =A cut . 

For a Gaussian density fluctuation field 8, weak lensing by 
large scale structure (e.g. Bartelmann & Schneider 2001; Miralda- 
Escude 1991; Blandford et al. 1991; Kaiser 1992; Jain & Seljak 
1992) produces a Gaussian magnification distribution. In fact, as 
long as 8A = A — 1 <C 1, the magnification of a source at redshift z 
can be approximated as 



A(z) = 1+SA(z)~1+2k(z) 



(54) 



i.e. to first order the magnification fluctuation is just twice the con- 
vergence K, which itself is a projection of the density contrast 8. 
The distribution of the magnification fluctuations 8A will then be 
Gaussian, with mean zero and a dispersion 0\(z) which depends 
on the source redshift and, albeit weakly, on cosmology. Typical 
values for Oa run from ~2x 10~ 3 for z = 0.05 to ~ 0.44 at z = 7.5 
(cf. Bartelmann & Schneider 2001). 

A convenient choice for A cut is A cut = 1 + 1.5<7 A (z), yielding 
Acut ~ 1.5 — 2 for the redshift range of interest in this paper. We 
model the probability distribution for A < A cut as 



p(A,z)=H(z) exp[-(A-A) 2 /2o A (z) 



(55) 



where the precise location of the peak A and the amplitude H(z) are 
determined by the normalization and flux conservation conditions 
[Eq. (^8|)] on the combined (weak plus strong lensing) probability 
distribution. 

The magnification probability for isothermal galaxy mod- 
els has been derived by Peacock (1982) and Vietri & Ostriker 
(1983), and, for more complicated galaxy models, by Blandford 
& Kochanek (1987) and Wallington & Narayan (1993). In this pa- 
per, we integrate Eq. ( j53| ) also for lenses with NFW density profile 
(see Sect. 0), and we describe the low-magnification distribution 
according to Eq. (p5|). Numerical results are presented in Section^ 

Let us now turn to the magnification bias on a flux-limited 
source sample. The integrated counts above a flux density threshold 
S v of sources with a comoving luminosity function 5>(L,z) can be 
written as (see e.g. De Zotti et al. 1996): 



N(S V 



( Z0 dz / dL<t>(L, 

JO JL min 



(1+z) 2 dz 



where the luminosity distance £>l(z, 1) is computed from Eq. ( p2| ) 
with (3=1 and 



W(V) = 431(1 +z)r z (z)S v 



L(v) 



L[(l+z)v] 



(57) 



Note that the comoving radial coordinate r [and its element dr in 
Eq. (po])] must not be confused with the dimensionless angular di- 
ameter distance r\ (z) [Eq. (p7|))]. 

The luminosity function modified by the magnification bias 
reads (e.g. Pei 1995): 



*'(M = f 

Ja„ 



p(A,z) L 
dA iy ' ; <j> -, 
A \A' 



(58) 



(56) 



Allowance for the effect of lensing on counts is made by replacing 
*'(L,z) to *(L,z) in Eq. M). 



6 SOURCE COUNTS IN THE SUBMILLIMETRE 
WAVEBAND 

In this paper, the "unlensed" galaxy counts in the submillime- 
tre waveband are taken from the model presented by Granato et 
al. (2000), which is in good agreement with the available SCUBA 
data at 850,um (Blain et al. 1999; Smail et al. 1999). 

In this model the rate of formation of spheroids at high redshift 
is estimated exploiting the (i) QSO Luminosity Function and (ii) 
observational evidence leading to the conclusion that high redshift 
QSOs did shine in the core of early type proto-galaxies during their 
main episode of star formation. In their scenario the star formation 
is more rapid in more massive objects, ranging from ~ 0.5 to ~ 2 
Gyr when going from more massive to less massive objects. This 
Anti-hierarchical Baryonic Collapse is expected to occur in Dark 
Matter halos, when the processes of cooling and heating is con- 
sidered. The larger the dark halo and the enclosed spheroid mass 
are, the shorter will the gas infall and cooling times be, resulting 
in a faster formation of the stars and of the central black hole. The 
star-formation process and the quasar shining phase proceed un- 
til powerful galactic winds are caused by the quasar itself, which 
occurs at a characteristic time when its luminosity becomes high 
enough. Also Monaco et al (2000), in order to account for the ob- 
served statistics of QSOs and elliptical galaxies in the framework of 
hierarchical structure formation, introduced a a time delay decreas- 
ing with mass between the beginning of the star formation and the 
QSO bright phase. 

The spectroscopic evolution of galaxies adopted here is based 
on the model GRASIL (Silva et al. 1998). GRASIL includes: (0 
chemical evolution; (ii) dust formation, assumed to follow the 
chemistry of the gas; (Hi) integrated spectra of simple stellar pop- 
ulations (SSP) with the appropriate chemical composition; (iv) re- 
alistic 3D distribution of stars, molecular clouds (in which stars 
form and subsequently escape) and diffuse dust; (v) radiative trans- 
fer computation in this clumpy ISM and dust temperature distribu- 
tion determined by the local radiation field. 

With these ingredients, the evolving luminosity functions (LF) 
at various wavelengths in the millimeter and submillimetre wave- 
bands are evaluated numerically, and they turn out to be signif- 
icantly different from models of Pure Luminosity Evolution for 
the 60,um LF of IRAS galaxies (Saunders et al. 1990), properly 
rescaled to the wavelengths of interest. The IRAS galaxy LF, which 
is based on an empirical model describing the evolution in a para- 
metric way, was used by Blain (1996) to obtain galaxy counts in the 
submillimetre wavebands, which were then used for estimating the 
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Figure 1. Magnification of an extended source as a function of the offset r 
between its center and the projection of the lens center in the source plane. 
Here the lens potential has no ellipticity. See text for details. 




Figure 2. Contour plot of the magnification for an extended source when 
g = 0.1. The projection of the lens center in the source plane is at the origin 
of axes. The levels are at A = 1.5,2,4, 6, 8, 10 from outside in. See text for 
details. 



incidence of gravitational lensing on source counts. However, the 
model by Granato et al. implies steeper source counts, nearly expo- 
nentially decreasing at bright fluxes, so that, as explained in section 
(Q), the effect of gravitational lensing is expected to be more im- 
portant. 

All source counts obtained in this paper, with and without 
magnification bias, the spheroids include elliptical galaxies as well 
as bulges of Sa galaxies, and we followed the formalism of Granato 
et al. (2000) on the source properties. 



7 MAXIMUM MAGNIFICATION FOR EXTENDED 
SOURCES 

As discussed by Peacock (1982), in the case of extended sources 
the magnification cannot be arbitrarily large. Correspondingly the 
magnification distribution must be cut-off at large A as: 



p(A) °c exp 



(59) 



where the cut-off magnification A max depends on the physical size 
of the source. 

Since, in some applications, e.g. for the estimate of the in- 
fluence of lensing on counts of sub-millimeter sources (Blain 
1996), the results are sensitive to the adopted value of A max , 
the approximated expression derived by Peacock (1982) [A max = 
70(DHq /c) (d /kpc) ~~ 1 , where d is the physical radius of the source, 
D is the angular diameter distance of the source and c/Hq is the 
Hubble radius] may not be sufficient. 

The morphology of strongly lensed sources indicates that most 
lenses are not circularly symmetric (e.g. Narayan & Bartelmann 
1997 and references therein). Therefore, to estimate the maximum 
possible magnification for a source of physical radius r at redshift 
z s , we considered in general an elliptical lensing potential due to 
a quasi isothermal sphere (Blandford & Kochanek 1987). Defining 
polar coordinates r and in the image plane, with origin at the 
center of mass, and measuring 6 from the major axis of the ellipse, 
the deflection potential may be written as: 



«P(r, 9) = 9 E J {s 2 + r 2 ) ( 1 - g cos 26) 



(60) 



Here g is the ellipticity parameter, s the core radius and 9g is de- 
fined by Eq. @. The results for an extended source are weakly 
dependent on s, within a relatively broad interval; therefore we sim- 
ply set s — in the following. With this potential we have computed 
the expected magnification of a pointlike source A(y), as a function 
of its position y in the source plane. This has been done by means 
of the ray-shooting method (e.g. Schneider et al. 1992, pag. 304). 
Then the magnification of an extended source with brightness pro- 
file I(y), as a function of the position ?£ of its center, is given by 



Ae{9e) = 



JI(y)A(y)d z y 
Jl(y)d 2 y 



(61) 



For the brightness profile, we have used either a de Vaucouleurs 
law 



log/(fl) = log/, - 3.33 [(R/R e ) [/4 - l] 
or the Hubble profile 
I{R)=I /(l+R/R„) 2 . 



(62) 



(63) 



Observed profiles of spheroidal galaxies are well reproduced by 
both functional forms over a wide range of R, provided that R e ~ 
11 R„ (Mihalas & Binney 1981). 

We find that, adopting this relationship between the scale- 
lengths, the estimated Ae are very similar in both cases. Therefore 
in the following we present only results for Eq. ([62]). We assume in 
particular/?, = 5 kpc, typical for a bright elliptical galaxy (e.g. Kent 
1985). For a SIS, the maximum possible magnification is achieved 
when the lens is local, i.e. when D^ s = D s in Eq. (pOp . We set 
a v = 300 km s , which corresponds to a mass of IO'^Mq within 
25 kpc adopting the isothermal relationship M(< R) = 2a 2 R/G. 
With our adopted cosmology, Eq. ( |2c| ) yields 9^ = 2.6". 

Figure hi shows the corresponding magnification for an ex- 
tended source as a function of the offset r between its center and 
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Figure 3. SIS model: comparison of the high-magnification tail amplitudes 
a(z) for a flat ACDM model with £1 A = 0.7 (solid line), and a flat CDM 
model (dotted line), both normalized to the cluster abundance. The ampli- 
tude is plotted vs. source redshift. The dot-dashed line refers to a COBE- 
normalized flat CDM universe. 



the projection of the lens center in the source plane. Here we have 
set g = into Eq. (|sc|), i.e. we adopt a spherical potential, and the 
source is at z s = 3, but very similar plots are obtained for z s in the 
range 1 4- 4. As can be seen, the magnification is maximized when 
the source and the lens are aligned and A max ~ 26. If the lens is 
placed instead at z/ =0.5, then b = 1.2 arcsec and A m ax — 

13. 

However, as already remarked, strongly lensed sources are 
usually explained in terms of potentials with non vanishing ellip- 
ticity. A typical value for g could be 0.1. In this case the symme- 
try of Ae around r = is lost and the maximum magnification, 
Amax — 12, significantly lower than in the spherical case, occurs 
when the lens and the source are offset by ~ 0.7" along the major 
axis of the ellipse. The results are detailed in Fig. ^. 

We conclude that reasonable values of A max for extended 
sources are in the range 10-30. The lower value is a relatively con- 
servative lower limit, easily exceeded for a wide range of values for 
the relevant quantities, while the latter is obtained only under rather 
special conditions. 



8 NUMERICAL RESULTS AND DISCUSSION 
8.1 Effects from cosmology 

To illustrate the effect of the cosmological model on the magnifi- 
cation probability distribution p(A,z), we plot in Fig. ^| the ampli- 
tude a(z) of the high-magnification tail for two different cosmolo- 
gies, viz. a flat CDM and a ACDM model with Q, 0A = 0.7 and 
flom =0.3. We have adopted SIS profiles for lenses and a smooth- 
ness parameter oc s = 0.9. The two models are normalized to repro- 
duce the local abundance of rich clusters: Og = 0.56f2 0| j*' 47 (Viana 
& Liddle 1999). 
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Figure 4. The Sheth & Tormen (1999) mass function for a flat CDM model 
(dotted line) and for a ACDM model (solid line), at redshifts z = 0. 1,2, 
respectively, from top to bottom. 

The optical depth for a beam of light from a source due to lens- 
ing is proportional to the number density of deflectors multiplied 
by the cross section for a given magnification, integrated along the 
line of sight. Two competing effects are therefore relevant. On one 
hand, the path length to a source is larger for a ACDM model. On 
the other hand, the structure formation histories within the stan- 
dard hierarchical clustering scenario are also different for the dif- 
ferent cosmologies. If we normalize the models to reproduce the 
local cluster abundance, the density of lower mass objects is lower 
in the ACDM model, because the mass function is flatter (e.g. Eke 
et al. 1996). This is shown in the upper panel of Fig. ^, where the 
mass function [Eq. ([j])] is plotted as a function of the halo mass at 
z = 0. The number of low-mass objects keeps increasing with red- 
shift above the ACDM model, while the opposite is true for high 
masses. 

The next question is then: which is the mass range contributing 
most to the optical depth? In Fig. ^ we show as a function of the 
lens mass, both the cross section for magnifications A > 2 and the 
lens mass function for given values of the source and lens redshifts. 
A SIS profile is adopted. It may be noted that the magnification 
cross section is similar for the two models (although slightly higher 
for the flat CDM model). The product of the two functions peaks 
at masses between 10 11 and 1O 12 M for both a flat CDM and a 
ACDM model. 

The effective mass of dark matter halos contributing most to 
strong lensing of a source located at z s can be estimated as: 

W -./ dZl / dMM ^#' (M) 

where P(A,z s ) is the probability to have magnification > A. The 
inner integral, d{M) /dzi, is plotted as a function of z\, for z s = 5, a 
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Figure 5. Number density of dark-matter halos (from the Sheth & Tormen 
mass function), and cross section for magnification A > 2 by a SIS lens, 
in arbitrary units, as a function of the lens mass in 1/ . . Dotted lines cor- 
respond to flat CDM, solid lines to ACDM, both normalized to the cluster 
abundance. 



Figure 6. Contributions from different redshifts to the effective lens mass 
(see text) for a fixed configuration of the lens system, SIS lenses, and a 
ACDM model. 



SIS lens profile and a ACDM model, in Fig. ^| As illustrated by this 
Figure, the maximum contribution to the magnification probability 
comes from the mass range (I0 ll -10 12 Mq) for which space den- 
sities implied by a CDM model are appreciably higher than in the 
case of a ACDM model, in the relevant redshift interval, if the mod- 
els have to be consistent with the observed cluster abundance. This 
over-compensates for the larger path length to a source in a ACDM 
model and explains why the probability distribution of strong mag- 
nifications has a lower amplitude for such model. 

It may be noted (Fig. ^ that the difference between the two 
cases decreases with increasing source redshift, since also the 
ACDM model approaches an Einstein-de Sitter universe at high 
redshift. 

Also shown in Fig. Q is the amplitude a(z) for a COBE- 
normalized standard-CDM model. The relatively large values for 
this quantity for any source redshift are unrealistic since they are 
due to a mass function of dark halos inconsistent with the cluster 
abundance. 



8.2 SIS versus NFW profiles 

Figure [7] compares the effect of SIS and NFW lens profiles on the 
magnification distributions for sources at z s = 4 and z s = 7, includ- 
ing the weak lensing effect for A < A cut . The magnification dis- 
tributions are obtained by integrating over lens masses the cross 
sections in the source plane, weighted by the mass distribution 
[Eq. (0)]. The weak-lensing regime, responsible for magnifications 
below A cut , gives rise to a Gaussian peak near A = 1 whose disper- 
sion increases with increasing source redshift. In the strong-lensing 
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Figure 7. Magnification distribution from a population of SIS (solid lines) 
and NFW lenses (dot-dashed lines) for sources at redshifts z s = 4 (lower 
curves) and z s = 1 (upper curves). 
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Figure 8. High-magnification tail of P(A). Magnifications are plotted on 
a logarithmic scale. The plots refer to populations of SIS (solid lines) and 
NFW lenses (dot-dashed lines) for sources at redshifts z s = 4 (lower curves) 
or at z s = 7 (upper curves). 



Figure 9. Magnification cross sections o"(A) (in square arcsec) for A > 2 
(upper curves) and A > 10, as a function of the halo mass. The sources are 
at z s = 5 and the lenses at z = 1 ■ The cross sections are plotted for SIS (solid 
lines) and NFW halos (dotted lines) 



regime, the probability distribution is obtained by solving Eq. ( |12| ) 
numerically and inserting the result in Eq. (^). The asymptotic be- 
haviour oc A~ 3 is reached earlier by the most distant sources. The 
high-magnification tail of Fig. ^ is shown in closer detail in Fig. |^ 
Note that the plotted distribution of magnifications has a discon- 
tinuity in A cut , i.e. at the transition between the weak and strong 
lensing regimes. 

The two density profiles lead to slightly different magnifi- 
cation distributions. In particular, the NFW lens is more efficient 
than the SIS for moderate magnifications (2 ^ A 4), and less ef- 
ficient for high magnifications. In fact, NFW lenses have smaller 
high-magnification cross sections than SIS lenses of equal mass, 
even if the average magnification is higher. This can be read off 
Fig. ^ where we plot the cross section for magnifications A > 2 and 
A > 10 as a function of the halo mass for the two models, keeping 
fixed the configuration of the system. For virtually all halo masses, 
o(10)nfw < <*(10)sis> while o(2) NFW > g(2)sis- Even though the 
latter relation fails for very small lens masses, it still holds when 
the cross sections are weighted with the appropriate mass function. 
As mentioned above, the bulk of the contribution to the magnifi- 
cation distribution comes from a limited mass range. The effective 
mass defined in Sect. ^7T] is nearly equal for SIS and NFW profiles, 
namely ~ 10 11_12 M e . Although (M) depends (albeit weakly) on 
the lens redshift, massive clusters never contribute substantially to 
the integrand in Eq. ( |53| ) because they are extremely rare. 

The SIS and NFW density profiles and the corresponding mass 
encompassed within the radius r for a 10 12 object are shown in 
Fig. hd. Virial radii are approximately equal in the two cases. For 
the NFW profile, the virial and scale radii are r2oo — 140 h kpc 



and r s ~ 15 A kpc, respectively, for h = 0.65; the concentration is 
c ~ 10. The con verg ence K for the two profiles is shown in the up- 
per panel of Fig. Illl as a function of halocentric distance. Unlike the 
surface density, the convergence depends also on the geometry of 
the lens system. Here, the lens is at z\ = 1 and the source at z s = 5. 
The two convergence profiles are quite similar, but the shear is also 
playing a fundamental role to determine the magnification distri- 
bution. In the lower panel of Fig. [Tl], we plot, for the two profiles, 
det~ J4 , i.e. the image magnification fj of an image with impact 
parameter r in the lens plane. Note that, even if some images are 
demagnified by lensing, the total magnification of the source is al- 
ways > 1 (e.g. Schneider et al. 1992). The critical curves behave 
differently. Even though the NFW profile has a singular core, it has 
tangential and radial critical curves (Bartelmann 1996), while the 
SIS has only a tangential critical curve (whose caustic degenerates 
to a point for all axially symmetric lenses). Consequently, the max- 
imum image number is two for SIS and three for NFW lenses. 

Figure [IT] shows that the tangential critical curve of a SIS oc- 
curs at a larger radius than both critical curves of an NFW halo with 
equal mass. This means that high magnifications are favored in the 
SIS model, because the corresponding cross sections in the source 
plane is larger. On the other hand, the SIS profile yields lower total 
magnifications when the source lies well outside the outer caustic 
in the source plane. The NFW cross sections for /j tot ~ 2 generally 
overcome the SIS ones, due to the fact that the lensing potential is 
less curved for the NFW than for the SIS profile, hence this also 
occurs for more massive halos. In Fig. [l^, we plot the equilibrium 
configurations for a 1O 14 M0 halo. The NFW profile has virial and 
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Figure 10. Upper panel: density profile for a dark matter halo of 
IO^/t'Mq at z = 1. The mass density is plotted versus the halocentric 
distance for SIS (solid line) and NFW halos (dashed line). Lower panel: 
mass encompassed by radius r. 
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Figure 11. Upper panel: convergence K as a function of halocentric distance 
for a halo with M^ a = 10 12 h 1 V/ . at z\ = 1. The sources are at z s = 5. The 
plot refers to a SIS (solid line) or a NFW halo (dashed line). Lower panel: 
Magnification as a function of impact parameter, involving convergence and 
shear of the gravitational field (see text). 



scale radii of r200 — 690/j kpc and r s ~ 100 h kpc, respectively, 
for h = 0.65 (the concentration is c ~ 6.5). 

In Fig. |l3| we plot the convergence and the magnification for 
a single image as a function of its location. After adding all image 
magnifications, inverting to find the cross section a(/j) vs. source 
position, and integrating over all lenses up to the source redshift, 
this explains the enhancement of the NFW probability for low mag- 
nifications with respect to the SIS model, shown in Fig. ^, and the 
opposite effect for large magnifications. However, we see that the 
probability distributions are quite similar and the effects of lens- 
ing on a source population turn out to be nearly identical for lenses 
with the two density profiles. 

We finish with a cautionary note. The simulations resulting 
in the density profile of Eq. ( pB] ) did not have sufficient resolution 
for halo masses < lO lo fc _1 M . We are therefore extrapolating the 
validity of this density profile as well as the relation between halo 
concentration and mass to a mass range where higher-resolution 
Af-body simulations would be required. Furthermore, as Porciani & 
Madau (2000) pointed out, in order to obtain agreement with the 
data on image separations of QSOs one should consider baryonic 
cooling in dark matter halos, able to transform NFW halos into 
isothermal distributions for masses smaller than some threshold. In 
such a picture, DM halos are modeled as NFW halos only above a 
threshold that is certainly well above the smallest mass that we are 
considering here. 



8.3 Effects on the source counts 

We can finally estimate the strongly lensed source counts, and com- 
pare them with the "unlensed" and weakly lensed counts, in a given 
wave band, based on the model by Granato et al. (2000). Using the 
"corrected" magnification distribution p(A) for extended sources, 
Eq. (^), we compute lensed counts at 850,um, using SIS and NFW 
lenses. In Fig. Il4j, the solid line shows the integral source counts 
that we expect at 850/7111 from the source distribution described in 
Sect. ^, ignoring lensing. The dot-dashed line includes only weak 
lensing, using the low-magnification tail of the magnification dis- 
tribution p(A,z)- Since the latter is modeled as a Gaussian with 
small dispersion around the mean A = 1, we can see from Fig. [l4| 
that weak lensing by large-scale structures has very little effect on 
the integral source counts, even though the variance of the distri- 
bution increases with source redshift. The weakly lensed counts 
are therefore quite similar to the unlensed ones given by Granato 
et al. (2000). Even the weakly lensed counts, however, fall above 
the unlensed counts where the number-count function falls most 
steeply. 

The effects of strong lensing from a Sheth & Tormen mass dis- 
tribution of SIS and NFW lenses is plotted in Fig. |l4| as short- and 
long-dashed lines, respectively; for both density profiles, following 
the results of section we show in Fig. [H] the results obtained 
plugging either A E;lnax = 10 (lower lines)or A E . max = 30 (upper 
lines). 

We can see that the contribution of strong lensing from a SIS 
model is here of the order 10~ 2 in the flat part of the counts, while 
it dramatically raises up in the steepest region, overcoming the 
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Figure 12. Upper panel: density profile for a dark matter halo of 
lO 14 /!-'M at z = 1. The mass density is plotted versus the halocentric 
distance for SIS (solid line) and NFW halos (dashed line). Lower panel: 
mass encompassed by radius r. 

weakly lensed counts. This is due to the very strong magnifica- 
tion bias for these sources. First, the source redshift, z — 5, is quite 
high. The probability for a source to undergo a lensing event with 
high magnification increases with increasing source redshift. Sec- 
ond, and importantly, the steep source counts, discussed in Sect. ^ 
and by Blain (1996), provide a huge reservoir of source to be mag- 
nified above the flux limit of the observation. 

At flux densities around lOOmJy, the counts will be domi- 
nated by sources magnified by a factor A ~ 10 — 20. As we have 
seen, such quite high magnifications are less probable for the NFW 
lenses than for the SIS population, which explains the relative val- 
ues of strongly lensed source counts from NFW and SIS lenses at 
flux densities £ lOOmJy of Fig. |l4[ At higher flux densities, lensed 
sources are depressed by the effect of having taken a maximum am- 
plification AE,max = 

30. 

The substantial magnification bias at flux densities around 
lOOmJy allows to detect sources otherwise not bright enough for 
detection. This result may have very interesting consequences on 
the expected number of SCUBA sources to be found in the whole 
sky maps produced by the incoming Planck Surveyor Satellite, op- 
erating at nine frequency channels between 0.3 and 10 mm (Man- 
dolesi et al. 1998). The detailed study of this issue will be the sub- 
ject of a forthcoming work (Perrotta et al. 2001). 
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